On a class of vector fields with discontinuity of divide-by-zero type 

and its applications 



R. Ghezzi* A.O. Remizov^" 

SISSA/ISAS, via Bonomea 265, 34136 Trieste, Italy 
^ ! 
O 

(N 

Let n be a domain in M", n > 2, with coordinates x = {xi, . . . , Xn)- Let / G C^{Q, M), s > 2, 
iy^' be such that the equation f{x) = defines a regular hypersurface F C fi, i.e., at all points 

Q ! X eT the condition V/(x) 7^ holds. We consider vector fields of the type 
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w{x) = r{x)v{x), (i.i: 



hm 


r\x)D^{x) 


= 


0, 


Vx, 


er, 


lim 

X^Xt 


nx)D^ix) = 


0, 




hm/ 


dxi 



where V G C*(fi,]R") is a vector field and r is a positive real number. 

The divergence of W , denoted by Dy^-, is infinite or undetermined on the hypersurface F, 
^ ■ but it is a C"'~^-smooth function on fi \ F. Assume that the field W satisfies the following 

CN ! main conditions 

p 

O ; lim r{x)D^{x) = 0, lim f+'::^(x) =0, Vx, G F : ^(x.) = 0, V^. 

For simplicity we write these assumptions in the form 
X; r'D^\^ = 0, (1.2) 

^' rD^{x,)=0, f+iJ(x,)=0, Vx. gF:V(x,) = 0, V^. (1.3) 

OXi 

Conditions (II. 2p . (I1.3P are fulfilled, for instance, if the vector field W is divergence-free, i.e., 
D^ = for all X G f^ \ F. 

The field W is C*-smooth on fi \ F, but at points of F formula (II. ip gives a discontinuity 
of divide- by-zero type. Due to their large number of applications (e.g., in mechanics with dry 
friction and control theory, see [7J), discontinuous vector fields (or, equivalently, differential 
equations with discontinuous righthand sides) have been widely studied. However to the 
authors' knowledge, fields of the given type have not been studied yet. 

Although this problem seems at first sight rather theoretical and not natural, it is mo- 
tivated by a large number of applications. Indeed, many variational problems in differen- 
tial geometry and calculus of variations are characterized by Lagrangian (or Hamiltonian) 
functions that are smooth at all points except for a regular hypersurface F. The vector field 
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corresponding to the Euler-Lagrange equations of such problems is divergence-free and takes 
the form (11.11) . The simplest example is the equation of geodesic lines on the cuspidal edge 
embedded in the Euclidean space or on the plane with the Klein metric, that is used in the 
model of the Lobachevsky plane. 

The aim of this paper is to establish some general facts about vector fields of the form 
(II. ip under assumptions (11.21) . (11.31) that allow to infer some properties on the vector fields 
V. Applications to several concrete problems, that are interesting from different points of 
view, are provided. 

The paper is orginized as follows. In section [2] we prove several simple theorems about 
vector (and direction) fields of the form (11.11) under assumptions (II. 2p . (II. 3p without any 
special hypothesis on V. In particular, these results show the key role of singular points of 
the field V in the applications. 

In section [3] we give a brief survey of the theory of normal forms at non-isolated singular 
points of smooth vector fields. We restrict to the case where the components of the vector 
field belong to the ideal generated by two of them in the ring of smooth functions. As far as 
we know, the first work devoted to the analysis of local normal forms for such fields is due to 
F. Takens [20] ■ Later, the problem was deeply investigated in finite smooth [12], C°°-smooth 
[T7] . and analytic categories pT| fF2\r\ 

In this survey, we deal only with finite and C°°-smooth classifications, which are simpler 
than the analytic one. Notice that almost all the facts in this section were previously known. 
We try to present the subject so that it is not obscured by technical details and, at the same 
time, is sufficiently precise. We hope that the informed reader will tolerate trivial aspects 
while the reader unfamiliar with this subject will understand the main ideas and find all the 
omitted proofs and technical details in the cited literature (see also [9l [11], [191 l2l] devoted 
to similar problems). 

In the last section we apply the results to the problem of geodesic flow generated by three 
different types of singular metrics on 2-surfaces. Firstly, we consider pseudo-Riemannian 
metrics, i.e., metrics that degenerate (change their signature) on a curve, see also [T5] . 
Secondly, we analyse metrics of Klein type, that are positive definite but have a singularity 
of divide-by- zero type, see [l^. Finally, we consider almost-Riemannian metrics, i.e., metrics 
whose orthonormal frames are pair of vector fields that are coUinear on a regular curve, see [1]. 
Two another examples can be found in [T5[ IT6] . 

2 Basic Theorems 

Integral curves of the fields W and V coincide at all points a; G fi \ F. At the same time 
the field V is more suitable for analysis, since it is smooth on the whole domain Q while the 
field W is discontinuous on the hypersurface F C fi. Our concern is to pass from the initial 
vector field W to the vector field V. 

— * 

Theorem 1 Condition (11.21) holds true if and only if F is an invariant hypersurface ofV. 
The function f is a first integral of V if and only if 

fD^ix) ^ D^{x). (2.1) 

Assume f to he a first integral ofV and let condition (11.31) holds true. Then D^(a;*) = for 
every x^ &V such that V{x^) = 0. 



^ Remark that the finite smooth classification is based on the general results by V.S. Samovol |18l . and the analytic 
classification is based on the general results by A.D. Bryuno H] — [6] and J.C. Yoccoz I23| . 



Proof. Using the formula of divergence in Cartesian coordinates, for every point in i7 \ F 
we get 

r+'D^ = fD^ + r^'L^iD = fD^ - rL^f, (2.2) 

where L^ denotes the Lie derivative along the vector field V. All terms in the right hand 
side of equation (12. 2p are C*~^-smooth on Q. Hence, taking the limit as x tends to x* G F 



it follows that f+^D 
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is equivalent to L-r^fl = 0. 



As concerns the second statement, for every point in fi \ F, we have 
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If identity (^^ holds on 11 \ F, then ([23D implies L^f\ 
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(2.3) 
0. By continuity it follows 



L^f = on r2, i.e., / is a first integral of V. Conversely, assume L^f = on 17. Then, by 
(12. 3p we get {f^D-^ — -Dy)|^,p = 0. Thus, by continuity {D-^ is continuous on f2), identity 

(12. ip on VL follows. Finally, combining the first equality in condition (II. Sp and identity (12. ip . 
we get the last statement of the theorem. q 

Corollary 1 Assume that condition (II. 2p holds. Let 7 he an integral curve of either V or 
W passing through the point x^ G F. IfV{x^) 7^ 0, then in a neighborhood of x^ the curve 7 
lies entirely in the hypersurface F. 

Theorem [1] and Corollary [1] explain why singular points of V play an important role. 
Indeed, in many applications it is necessary to find integral curves that intersect the invariant 
hypersurface F but do not belong entirely to F. Hence such integral curves intersect F only 
at singular points. The next theorem establishes a relation between the eigenvalues of the 
linearization of the vector field \^ at a singular point x^, eT. As we shall see, in many cases 
such a relation is resonance. 




Figure 1: Three examples of phase portraits of the vector field V. a) The case V{x^,) ^ 0, and all integral 
curves belong to F. b), c) The case V{x^,) = 0, and all integral curves (except only one) do not belong to F. 



Theorem 2 Let x^ E T be a singular point of the field V and Ai, . . . , A„ be the eigenvalues 

of the linearization of V at x^,. If conditions (ll.2p and (I1.3P hold true then there exists 
j G {1, . . . , n} such that 

Ai + • • • + A„ = rXj. (2.4) 

There exists an eigenvector corresponding to Xj which is transversal to F at x^. The spectrum 
of the linearization of the restriction V\„ at x* is {Ai, . . . , A„} \ Xj. 
If f is a first integral of V then Xj = 0. 



Proof. Differentiating the identity f^^^D^ = fD^ — rL^f (see (12.21) ) with respect to Xj 
we get 

[r + l)fD^^— + f+ -—^ = £).-_ + /_-J^ - r( -— , Vf ) -r( V, V—- 

OXi OXi OXi OXi \OXi I \ OXi 

where the triangle brackets denote the standard scalar product of vectors. 

The last equality holds for all x G f2 \ F and its right hand side is C"*~^-smooth on f2. 
Taking the limit as x tends to x* G F such that V{x^ = and using (II. 3p . we get 
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0, i = 1 n. 



Since r 7^ 0, the last system can be written in matrix form as Ae = pe, where A = (|j^) I 

is the matrix of the linearization of V at the singular point x*, the vector e = V/(x^,), and 
the number p = r~^D^{x^:). By hypothesis, V/(x*) 7^ 0, hence p is an eigenvalue of the 
linearization of V at x* with corresponding eigenvector V/(x,,). Let j G {1, . . .n} be such 
that p = Xj. Notice that the divergence of a vector field at any singular point coincides with 
the trace of the linearization of this field at that point. Thus Xj = r^^(Ai + • ■ ■ + A„) which 
leads to equality (12. 4p . Clearly, e is transversal to F at x,,,, whence Xj does not belong to the 
spectrum of the linearization of the restriction V\ . 

To prove the last statement recall that if / is a first integral of V, by Theorem [T]D^(x*) = 
follows. Since D^(x*) = Ai + ■ ■ ■ + An, we have equality (12. 4p with Aj = 0. q 

Remark 1 TheoremsUl andl^ hold true not only for r > 0, but also for r < 0. 

Let us illustrate two examples in M^ with coordinates {x,y, z). 
Example 1. Let V be the vector field 

X = X, y = y, z = z. 

Then D^{x, y,z) = 3 and the unique singular point of V is the origin. The spectrum of the 

linearization at the origin is (1, 1, 1). The field V has no non-constant first integrals, but it 
has the family of integral planes ax + by + cz = passing through the origin. 

Consider the vector field W given by formula (I LI I) with /(x, y, z) = ax + by + cz. Then 
^w ~ (2 ^''^)f~^ ^iid f^D^ = (3 — r). Condition (IL2p is satisfied, but (ILSp is fulfilled only 
if r = 3. In the case r = 3, we have relation (12.40 with any index j = 1, 2, 3. 

Example 2. Let V be the vector field 

X = 2x, y = y, z = 0. 

Then -D^(x, y,z) = 3 and the set of singular points of V is the z-axis. The spectrum of the 
linearization at any singular point is (2, 1, 0). The coordinate function 2; is a first integral of 
V and there is a family of integral surfaces given by x — cy"^ = 0, as c varies in M. 

Consider the vector field W given by formula (ILip with /(x, y, z) = z. Then f^D^ = 3 
and condition (II. 3p is not satisfied. This corresponds to the last claim of Theorem [H Indeed, 
since / is a first integral of V, condition (IL3p would imply -D^(0, 0, z) = 0. 

Let now f{x,y, z) = x — cy"^. Defining W as in (II. ip . we get f^D^ = (3 — 2r). Hence 
condition (II. 3p holds true if and only if r = 3/2 and the relation (12. 4p is valid with Aj = 2. 



Finally, consider the vector field W with f{x,y, z) = y. Then fD^^ = (3 — r), condition 
fll.Sp holds in the case r = 3, and we have (12.41) with Xj = 1. 

Sometimes it is more natural to consider direction fields rather than vector fields. Recall 
that given a vector field V, the direction field x associated to V is the class of vector fields 
(fV, where (p G C'^{Q) never vanishes. Theorems [H [2] are valid for direction fields. 

Theorem 3 Let ip G C^{Q) and ^plx) ^ for every x E Q. Then TheoremsUl hold true 
if in equation fll.ip we replace V by ipV . 

Proof. It is necessary and sufficient to prove that the main assumptions (II. 2p and (11.30 
are invariant with respect to multiplication of the vector fields V (and consequently W) by 
a C"^-smooth scalar factor (p y^ 0. Indeed, D -^^ = pD^ + f'^'L^p. Hence we get 
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These expressions show that conditions ( II. 2p . ( II. 3p hold true for the vector fields pV , ipW. 

D 

3 Fields with non-isolated singular points 

We are interested in studying vector fields V of the form 

C.=v, r] = w, (i = aiV + /3iW, i = l,...,l, (3.1) 

where ctj, /3i and v, w are C°°-smooth functions of the variables ^, ri,(i, ■ ■ ■ , Ci- Such a kind of 
fields occurs in many problems, for instance, in studying implicit differential equations (see 
next example) and slow-fast systems. 

Example 3. Consider the family of first-order implicit differential equations 

dx 
F{t,x,p)=e, P=^, (3.2) 

depending on the real parameter e not necessarily small. One effective approach (which 
goes back to Poincare) consists of lifting the multi- valued direction field defined by equation 
( 13. 2 p on the (t, x)-plane to a single-valued direction field x defined by equation (13.21) in the 
{t,x,p)-spa.ce. 

Geometrically, x is the intersection between the contact planes dx = pdt and the planes 
tangent to the surfaces {F = e} with various e. This gives the Pfaffian system 

Ft dt + F^dx + Fp dp = 0, pdt — dx = 0. 

Whence the direction field x corresponds to the vector field V given by the formula 

i = Fp, x = pFp, p=~{Ft+pF,), (3.3) 

where a dot over a symbol means differentiation with respect to the independent variable 
playing the role of time. The field ( 13. 3 p has the form (13. ip with / = 1, where ^ = t, rj = p, 

C = X and V = Fp, w = —{Ft + pF^). 



In this section we recall smooth local normal forms of fields f l3.ip at singular points. 

The components of V belong to the ideal I = {v, w) generated by two of them in the ring 
of germs of C°°-smooth functions (this property is invariant with respect to the action of 
diffeomorphisms of the phase space). The set of singular points of V is defined by the two 
equations v = w = 0. The spectrum of the linearization of V at any singular point contains 
at least / zero eigenvalues, i.e., it is (Ai, A2, 0, . . . , 0). 

Consider the germ of \^ at a given singular point. Without loss of generality we may 
assume the singular point to be the origin of the phase space. From now on, we will always 
assume that ReAi,2(0) 7^ 0, whence the set of singular points of V is the regular center 
manifold of dimension /, denoted by W^. The eigenvectors with zero eigenvalue are tangent 
to W^ and the eigenvectors (if they exist) corresponding to Ai,2(0) are tangent to the plane 
dCi = aid$, + iSidT], i = 1,. . .,1. 

It is convenient to choose local coordinates (^, ri,(i, ■ ■ ■ , (i) such that W^ ={^ = ''7 = 0} 
and the linear part of V at is in normal Jordan form. Then there exist C°°-smooth 
functions f 1,2 and Wi^2 such that v = ^f 1 + riV2 and w = ^Wi + rjW2. The eigenvalues Ai,2 
at various singular points continuously depend on the variable C, = (Ci, • • • ,0); which is a 
local coordinate on W^ (this dependence is C°°-smooth at the points where Ai 7^ A2). In the 
following we will always work using such coordinates. 

3.1 Normal forms: the non-resonant case 

We shall say that k functions U^^\^, r],Q, i = 1, . . . ,k, are independent by ( at the point if 
their gradients with respect to the variable ( = (Ci, . . . , 0) at are linearly independent. It is 
not hard to see that if a function f/ is a first integral of V then its partial derivatives U^ and 
Urj vanish at 0. Hence the number of first integrals of V independent at is not greater than 
/. On the other hand, the restriction of V to W^ is identically zero, hence by Shoshitaishvili's 
reduction theorem |2], the germ of \^ at is orbitally topologically equivalent to 

i = ^, V = ±V, = 0, i = l,...,l. 

The trivial equations Q = suggest the existence of / independent first integrals of V. If 
a /-uple of smooth first integrals U^^^ , • • • , U^''^ independent by C at exists, the change of 
coordinates Q ^-^■ f/'-*-' brings the field V to the form 

^ = v, r] = w, = 0, z = 1, . . . , /, 

where v and w are smooth functions obtained from v and w by the above change of coordi- 
nates. 

The existence of / independent smooth first integrals is connected with the following 
concept. 

Definition 1 The relations 

PiXi + P2>^2 = 0, pi,2 e Z+, Pl+P2>l, (3.4) 

are called resonances of first type. The minimal number pi+p2 (i.e., pi and p2 are relatively 
prime) is called the order of the resonance (13.41) . 

Consider the germ of a smooth function U{^, 77, () at the point and its Taylor series with 
respect to the variables ^,77, i.e.. 



The germ of U is called N-flat (A^ G N or oo) by ^, tj if Up^p^{Q = for all Pi + P2 < ^• 

Lemma 1 If between the eigenvalues Ai^2(0) there are no resonances (13. 4 p up to order N G N 
inclusive, then there exist C^ -smooth functions U^^\ . . . , U^^^ independent by ( atO such that 
L^U^^\ ..., L^U^'-^ are N-fiat by^,r] atO. 

Proof. We prove the lemma assuming Ai,2(0) to be real and the Jordan form of the 
linearization of V^ at to be diagonal (the cases of complex eigenvalues or Jordan form with 
a second-order cell to be considered similarly). Then the germ of V has the form 

e = e(Ai(o)+^i(c) + ---) + ^(^2(0 + ■■■), 

V = e(^i(C) + ■■■) + ^(A2(0) + ^2(0 + ■■■), (3.6) 

Q = ah?ic) + ---)+vih^2\c) + ---), ^ = 1,...,/, 

where all functions tIi,2(C)? 'S'i,2(C)? ^1 2(C) are C°°-smooth and vanish at 0, and the omitted 
terms are C°°-smooth functions containing the factor ^ or rj. 

The idea is to look for the functions U^^\ . . . , U^'''' in the set of polynomials in ^, 77 with 
coefficients smoothly depending on (. Namely, consider a function U in the form (13. 5p with 
finite sum < pi + p2 < N and unknown coefficients MpipalO- Substituting this expression 
into L-^U and using (13. 6p . we get 

N , 

Pl+P2=0 ^ 

+ up,p,pie'-V'-'\v2 + ■■■)+ up,p,P2e'^'v'"^\w^ + •■■)+ 

Z^{e'^V'{h['^ + ■■■)+ e'rf^^\hi'^ + ■■■)))■ 



/ 



k= 



For L^U to be 1-flat by C,, rj we set the coefficients of the monomials ^ and rj equal to 
zero, that is, 

(Ai(0) + v,)u,, + w,u,, + V /.f )^ = 0, 

'7' (3.7) 



V2U,o + MO) + W2)Uoi + V/^f ^ = 0. 



fc=l 

Since the determinant di{() of the linear system (13. 7p with respect to the unknown variables 
tiio and Moi is a smooth function and di{0) = Ai(0)A2(0) 7^ 0, the solutions Uiq and Mqi are 
smooth in a neighborhood of 0. Notice that the functions 'Uio(C) and tioi(C) depend on the 
derivatives ^^, where Moo(C) is any arbitrary smooth function. 

Given n G {2, . . . , A^} consider the coefficients of the monomials ^P^r/^^, where pi +p2 = n. 
In order L^U to be n-flat, we set these coefficients to be identically zero, i.e., we solve the 
system 

(Pl(Ai(0)+t^i)+P2(A2(0)+M72))Mpip2 + (Pl + l)^2Mpi+lp2~l + (P2 + l)WlMpi-lp2+l = </'piP2> (3-8) 

where v^p^pa are polynomials of the coefficients Uo,p{Q and their first-order derivatives with 
a + (3 < n (with Uap{C) = if a < or /3 < 0). The determinant (in(C) of the linear system 



fl3.8p with respect to variables Up^p^ , pi + p2 = n, has the form 

dn{0= n (PlAl(0)+P2A2(0) + 5(C)), 

Pi+P2=n 

where 6{() is a smooth function vanishing at C = 0. The absence of resonances (13.4p up 
to order A^ imphes that dn{0) 7^ 0, whence in a neighborhood of the coefficients MpipalC); 
Pi + P2 = n, smoothly depend on the functions Ua/siC) and their first-order derivatives with 
a + /3 < n. 

Finally, let Mqo (C); • • • > "^oolC) be C°°-smooth functions independent at 0. For every index 
2 = 1, . . . , / we define f/*^*-* by solving systems (13. 7p and (13. 8p for n = 2, . . . , N with the initial 
function uqq = Uqq. By construction, U^^\ . . . , ?7*^'^ satisfy the required conditions. q 

Corollary 2 If between the eigenvalues Ai^2(0) there are no resonances (13. 4p up to order N 
inclusive, then the last I components aiV + (3iW of the field (13. ip can be assumed to be N-flat 
by ^,r] atO. 

Corollary |2] allows to get normal forms in C'^-smooth and C°°-smooth categories. As for 
the C'^-smooth category, in [18] the author defines the number 



N{k) 



^2^^^^max|ReA,2| 



+ 2, keN, 



min I Re Ai^2 
the square brackets denoting the integer part of a number. 

Theorem 4 // between the eigenvalues Ai^2(0) there are no resonances (13. 4p up to order 
N{k), then the germ of (13. ip is C'' -smoothly equivalent to 

^ = v, r] = w, Ci = 0, z = l,...,/, (3.9) 

where v,w are some new functions of ^,ri,(. If between the eigenvalues Ai,2(C) there are 
no resonances (13. 4p of any order for all ( sufficiently close to 0, then the germ of ( 13. ip is 
C°° -smoothly equivalent to (13. 9p . 

The proof of Theorem H] in the finite-smooth category is based on Lemma [Hand on general 
results from [TH] . The proof in the C°°-smooth category requires more advanced techniques 
(see [17] or [9]). Notice that if Re Ai^2(0) have the same sign the absence of resonances (13. 4p 
between Ai,2(0) implies the absence of resonances (13. 4p between Ai,2(C) for ciH C sufficiently 
close to 0. This is no longer true if Re Ai,2(0) have different signs, except for the special case 
when the ratio A = A1/A2 is constant on W^, i.e., at all singular points. 

As we shall see in the following, the normal form (13. 9 p can be further simplified. 

Definition 2 The relations 

PiAi+P2A2 = Aj, Pi,2eZ+, pi+P2>2, JG{1,2}, (3.10) 

are called resonances of second type. The number pi + p2 is the order of resonance. 

Clearly, a resonance (13. 4p of order n implies a resonance (I3.10p of order n -\- 1. In this 
section we assume the absence of resonances (13. 4p up to order A^ G N or 00. Hence a 
resonance (I3.10p of order < A^ holds if and only if the ratio A(0) = Ai(0)/A2(0) or its inverse 
belongs to {2, ... , A^}. Combining the results from [S], [12], [IZ], [IB], one gets the following 
theorems. 



Theorem 5 Let A; G N and assume that between Ai^2(0) there are no resonances ( 13.1 Op of 
order N{k) inclusive. Then the germ ofV at is C^ -smoothly equivalent to 

e = ai(C)e + «2(CK ^ = «3(c)e + a4(CK 6 = 0, 2 = 1,...,/. (3.11) 

Moreover, if Ai,2(0) are real and A(0) 7^ 1, the germ of V at is C'' -smoothly orhitally 
equivalent to 

e = A(C)e, n = ri, 6 = 0, 2 = 1,...,/. (3.12) 

Both statements hold true with k = 00 if between Ai,2(C) there are no resonances fl3.10p of 
any order for all ( sufficiently close to 0. 

— * 

Theorem 6 Assume that A(0) = n is natural. Then the germ of V at is C^-smoothly 
orbitally equivalent to 

e = A(C)e + V^(C)^", V = V, 6 = 0, t = l,...,L (3.13) 

If ip{0) 7^ 0, then ip{Q simplifies to 1; if (p{() has a zero of finite order s at the origin then 
f{C) simplifies to (^ . 

The normal forms (I3.1ip -( l3.13p show that in a smaU neighborhood of the phase portrait 
of V is rather simple and V has a smooth 2-dimensional invariant foliation given by the 
equation ( = c in normal coordinates. The restriction of V to each leaf C = c is a planar 
vector field with non-degenerate singular point: node, saddle, or focus. 

3.2 Normal forms: the resonant case 

Consider the case where between Ai^2(0) there is a resonance (13. 4p . i.e., there exist n,m eN 
relatively prime such that 

mAi + nA2 = 0. (3.14) 

In this case, the proof of Lemma [1] for N > n -\- m fails, since the determinant dn+m{0 of 
the linear system ( 13. Sp with pi -\- p2 = n -\- m vanishes at C = 0, and vector field ( 13. ip with 
resonance (I3.14p at may not have a Z-uple of smooth first integrals independent at 0. 

A simple illustration (with I = 1 and n = m = 1) comes from Example 3. Indeed, let 
be a singular point of the vector field V given by (13. 3p . Clearly, F is a first integral 
of V, and the derivatives Fp and Ft vanish at 0. Assume that Ai(0) + A2(0) = 0. Since 
Ai + A2 = -Dy = —Fx, we have ^^^(0) = 0, i.e., F is not regular at 0. Let F be another 
first integral of V. Then the integral curves of V are 1-graphs of solutions of the implicit 
equation F{x,y,p) = e with various e. Hence, the previous argument with F replacing F 
leads to the same conclusion. Thus the germ of V^ at admits no regular first integrals. 

The resonance (I3.14p generates two infinite sequences of resonances (I3.10p , namely, 

{mj + l)Ai + njX2 = Ai, mjXi + {nj + 1)A2 = A2, j G N. 

This suggests that the formal normal form contains infinite number of terms p^, ^p^ , rjp^ , 
where p = ^"^rj"' is called resonance monomial corresponding to (I3.14p . The central step in 
the derivation of normal forms in the resonant case is the following. 

Lemma 2 For any A; G N, the germ ofV at is C^-smoothly equivalent to 

e = ^(Ai(0) + $i(p,C)), 77 = r/(A2(0) + $2(p,C)), 6 = P^^(P,C), ^ = 1,...,/, (3.15) 



where $i2(p, C) '^^^^ ^i(P)C) '^'"^ polynomials in p = ^'"r]" 0/ degrees N{k) and N{k) — 1, 
respectively, with coefficients smoothly depending on (. 

Assume that \l/i(0,0) 7^ 0. Then for any ui, . . . ,ui G M the germ (13.151) has a smooth 
first integral U{p, Q such that 

$(p, C)f/p + ^lip, CWc, + ■ ■ ■ + v&Kp, C)f/o = (3.16) 

UpiO,0)=uu U^,i0,0)=U2,...,U^M0)=^h (3.17) 

where $(p,C) = m<l>i(p,C) +r2$2(p,C)- 

The proof of Lemma [2] is based on the general results in [18] and can be found in |,12| . 

From now on, we will always assume \l'i(0,0) 7^ 0. This hypothesis implies the existence 
of / — 1 independent first integrals f/*^^\ . . . , f/*-'-* given by solutions of fl3.16p with initial 
conditions fl3.17p corresponding to linearly independent {I — l)-uples {u2, ■ ■ ■ ,uJi). Applying 
the change of coordinates Q 1— )■ U^^\ i = 2,...,/, the vector field (I3.15P is C°°-smoothly 
equivalent to 

e = e(Ai(0) + $i(p,C)), v = v{H0) + Mp,O), Ci = P^i(p,C), 6 = 0, t = 2,...,l, 

(3.18) 
where <l>i,2(p, C) ^-^d \l/i(p, C) are smooth functions of p and ( (not necessarily polynomials 
in p like in IKT5\f ). $1,2(0, 0) = 0, and ^i(0, 0) 7^ 0. 

The first integral U{p, () given by the solution of f l3.16p with initial conditions ui = 1,ijJ2 = 
■ ■ ■ = ui = allows to simplify the form f l3.18p . Considering the restriction $(p, C)|pi/c = 
$(0,C), we analyse two cases: $(^^(0,0) 7^ 0, which is generic, or $(0,C) = 0, which occurs 
in the analysis of some concrete problems (for instance, when n = m = 1, this condition 
corresponds to divergence- free fields). 

In the first case, there exists a C°°-smooth change of coordinates that preserves the form 
f l3.18p and brings the first integral satisfying ( 13.1 6p with initial conditions Ui = 1,U2 = ■ ■ ■ = 
Ui = to the form U{p, C) = P + C^- Even if the form fl3.18p cannot be further simplified, 
the phase portrait of V can be described using the invariant foliation p + (^ = c, see [12] . 

Similarly, in the second case there exists a C°°-smooth change of coordinates that pre- 
serves the form f l3.18p and brings the first integral satisfying fl3.16p with initial conditions 
Ui = 1,U2 = ■■■ = uJi = to the form U{p, () = p. Using this fact, the normal form f l3.18p 
simplifies as follows. 

Theorem 7 If conditions \&i(0,0) 7^ and $(0,C) = in ( 13.1 5p hold, then the germ of V 
at is C°° -smoothly orbitally equivalent to 

^ = <, r] = -mr], (i = P, 6 = 0, 2 = 2, ...,/. (3.19) 

The normal form ( I3.19P with any n,m & N was established in the C'^-smooth category 
for arbitrary fc G N in [12]. It was previously proved by R. Roussarie for the partial case 
77, = m = 1 in C°°-smooth category [T7]. The techniques developed in [17] can be applied to 
establish the normal form fl3.19p with any n, m G N in the C°°-smooth category. However, 
to the authors' knowledge, this result is not published. 

Remark 2 Theorem [^ is not valid in the analytic case: the analytic normal form is ob- 
tained from the smooth normal form fl3.19p by adding some module, see I21\ \2^ . From the 
viewpoint of the general theory developed by A. D. Bryuno ^-/Sj'; ^^^■^ ^'^'^ be explained in 
the following way. The condition A for the formal normal form f l3.19p does not hold, since 
the third equation in fl3.19p has the form Ci = P instead 0/ Ci = 0. Moreover, the pair 
of nonzero eigenvalues Ai,2(0) lies in the Siegel domain, where formal normalizing series 
generally diverge. 
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Remark 3 The condition \l/i(0,0) ^ in LemmalE can be replaced by \l/j(0,0) 7^ for 
some i G {1,...,/}. This condition holds true for germs (13.1 p with resonance fl3.14p having 
generic {n + m)-jet. Moreover, in order to check this condition it is sufficient to bring only 
the {n + m)-jet of (13. ip to the form (I3.15p . 

The following example shows that the condition ^i(0, 0) 7^ is essential. 
Example 4. Consider the vector fields 

e = e, v = -v, c = o, (3.20) 

e = e, v = -v{i + ^v), c = e^c, (3.21) 

both having at each singular point the resonance (I3.14p with n = m = 1, whence $(0, Q = 0. 
Clearly, the plane {( = 0} is invariant for both the vector fields and it is transversal to the 
center manifold W^ = {C, = rj = 0} at the origin. If the germ of either (13.200 or (13.211) 
were C'^-smoothly {k > 2) orbitally equivalent to the normal form (13.190 . then (I3.19P had a 
C'^-smooth invariant surface transversal to the ^-axis, i.e., of the form ( = f{^,ri). On the 
other hand, substituting the Taylor expression (of the second degree) of the function /(^, rj) 
into equation C,f^ — rjf^ — ^rj = 0, it is not hard to see that (I3.19P can not have an invariant 
surface of the form ( = /(^, r/). 

Remark 4 If n + m is rather large and the ratio n/m is sufficiently close to 1, the inequality 
n+m > 2[(2A;+1) max {n/m, m/n}]+2 has solutions k E N. According to Theorem\^ for any 
such k the germ of (I3.19P is C^-smoothly orbitally equivalent to (I3.12p with \{C) = —n/m 
or, equivalently, to the field 

i = ni, r] = -mr], 4 = 0, i = l,...,l. 

4 Applications: geodesic flows on surfaces with singular metrics 

We start with some general consequences of the results in two previous sections and then 
apply them to several concrete problems connected with singularities of divide-by-zero type. 

Let W he a vector field of the type in (11.11) . where r 7^ 0, 1 and the smootlu vector field 
V has the form (13. ip . Assume that conditions (II. 2p and (11.30 hold true. Let be a singular 
point of V such that the linearization of V^ at has at least one non-zero real eigenvalue, 
i.e., the spectrum is (Ai, A2, 0, . . . , 0), where Ai G M \ {0}. 

By Theorem [2], we have equality (12.41) . which in this case reads Ai + A2 = rXj, where 
j = 1, 2, or Ai + A2 = 0. Each of these equalities defines the spectrum of V^ up to a common 
factor 0", i.e., it uniquely defines the spectrum of the corresponding direction field. In both 
cases Ai^2 G I^\ {0}, hence in a neighborhood of the set of singular points of V is the center 
manifold W^, codimiy^ = 2. 

Theorem 8 Assume W^ C T , then in a neighborhood of the following statements hold. 

1) There exists a smooth regular function g : T -^ M. such that W^ = {g = ^} o-nd 
V\^ = gVi\^, where Vi\^ is a smooth non-vanishing field on T. 

2) The spectrum of the linearization of V at any singular point is o"(l, r — 1, 0, . . . , 0), 
where a is a smooth non-vanishing function on W^. 

3) The field V is smoothly orbitally equivalent to one of the following normal forms: 



For simplicity, we always assume that smooth means C°°-smooth. 
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fl3.12p with \{z) = r -I if r > I and r - 1, (r - 1)"^ ^ N or r < I and r ^ Q, 

f l3.13p with X{z) = n if r — 1 or (r — 1)^^ is equal to n eN, 

f l3.19p if r — 1 = ~n/m, where n,m E N, and "^{{0, 0) 7^ for at least one index 
i = 1, ... ,1 in the preliminary form (13.151) . 

Proof. For the first statement choose local coordinates (^, rj, () such that the invariant 
hypersurface F is the hyperplane {^ = 0} and the center manifold W^ is the subspace 
{C, = 1] = 0}. Then the field V has the form 

where v, Wi^2 and a^, /3j are smooth functions of C,, rj, (, and Ai = f (0), A2 = w{0). Substitut- 
ing .^ = 0, the field V\^ has the form 1) = r]W2, (i = Vl^i''^2, i = I, ■ ■ ■ J- Setting the function 
g = r]W2, we get T D {g = 0} = W and V\^ = gVi\^, where the field Vi\^ has the form 

^ = 1, 6 = A, « = l,•••,^• 
As for the second statement, according to previous reasonings, at any singular point in a 
neighborhood of we have equality Ai + A2 = rXj, where j = 1, 2, or Ai + A2 = 0. From the 
hypothesis W^ C F it follows equality f l2.4p with Xj = is impossible. Indeed, by Theorem [2] 
the spectrum of the linearization of the restriction V|j, at is (Ai, A2, 0, . . . , 0), where the 

number of zero eigenvalues is less by 1 than in the spectrum of V, i.e., is equal to / — 1. On 
the other hand, the inclusion W^ C F implies that the spectrum of the linearization of the 
restriction V\^ contains / zeros. Hence we have equahty Ai + A2 = rXj, with j = 1 or 2. 
Without loss of generality one can put j = 1, then A2 = (r — l)Ai. Since the last equality 
holds identically at all points in W^, the spectrum is (T{l,r — 1,0,..., 0) with a smooth 
non- vanishing function a. The third statement follows from Theorems [5]-[7] and Remark [1] 

D 

Each of the applications in this section will cast in the following situation. 
Consider the Euler-Lagrange equation 

|l,-L. = 0, p=|, (4.1) 

with Lagrangian L(t,x,p), where t,x eM.. In the (t, a;,p)-space equation (14. ip generates the 
direction field x corresponding to the vector field W given by 

i = Lpp, X = pLpp, p = Lx — Ltp — pLxp, (4-2) 

where the dot over a symbol means differentiation with respect to an independent variable 
playing the role of time. 

Lemma 3 At all points of the {t,x,p) -space where L is smooth the identity D^ = holds. 
Consequently, at all singular points of the vector field W where L is smooth, the spectrum 
of the linearization ofW has resonance Ai + A2 = 0. The same statements are valid for the 
corresponding direction field x- 

Proof. The identity D^ = is due to simple calculation. The field (14. 2 p belongs to 
the class of vector fields of type (13. ip , where the generators of the ideal / are v = Lpp and 
w = Lx — Ltp — pLxp. Hence the spectrum of the linearization of W at any singular point is 
(Ai, A2, 0). The equality Ai + A2 = for the field W follows from the equality D^^ = 0. The 
same equality for the fields (pW follows from the identity D ■^J^ = (pD^ + Ly^ip. q 
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In the applications below we deal with the case when the Lagrangian is smooth at all 
points of the (t, x, p)-space except for the the points of some regular surface F = {/ = 0} and 
the components of the field W given by formula (14.21) are fractions with common denominator 
/'', r > 0. Thus the field W is connected with some smooth field V by the formula (11. ip . 
From the identity D^ = (Lemma [3]) it follows that conditions (II. 2p and (II. 3p will be 
always satisfied, hence Theorems [T]-[3] and [8] are valid. 

4.1 Pseudo-Riemannian metric 

Consider a surface S with a system of coordinates (t, x) and a pseudo-Riemannian metric 

Q{dt, dx) = a(t, x) dx^ + 26(t, x) dxdt + c(t, x) dt^ (4-3) 

with smooth coefficients a, 6, c. The quadratic form Q is positive definite on an open domain 
S G S (which is called elliptic), indefinite on some other open domain H G S (which is 
called hyperbolic) , and degenerate on the curve A = {A = 0}, where A = b"^ — ac is the 
discriminant of the form Q. The curve A separates the domains S and Ti, every point of A 
is said parabolic. 

Example 5. Let 5* be a smooth surface embedded in the 3-dimensional Minkowski 
space, i.e., the 3D affine space with Cartesian coordinates {x, y, z) endowed with the pseudo- 
Euclidean metric ds"^ = dx"^ + dy"^ — dz^. A pseudo-Riemannian metric is induced on S by 
the metric ds^ in the ambient space. Denote by Cp the light cone in the 3D tangent space 
at the point P = {x, y, z) given by the equation dx"^ + dy"^ — dz"^ = 0. Then three possibilities 
arise: either the tangent plane to S" at P does not intersect Cp (then P G £^), or it intersects 
Cp along a pair of lines (then P G "H), or finally it intersects Cp along a unique line (then 
P is parabolic). 

For instance, if S" is a Euclidean sphere (x^ + y'^ + z'^ = r^), the parabolic points form two 
circles z = ±r/\f2, which separate S into two elliptic domains {£ : \z\ > r/\/2) and one 
hyperbolic domain ("H : \z\ < r/^/2). Geodesies on Euclidean spheres and ellipsoids in 3D 
Minkowski space are well-studied, see e.g. [8], [10]. 

Consider geodesies generated by the pseudo-Riemannian metric (14. 3 p in a neighbourhood 
of a parabolic point. Their 1-graphs are extremals of equation (14. ip with L = a/P, where 
P = a{t, x)p'^ + 2b{t, x)p + c{t, x). Then the vector field W given by formula (14. 2 p reads 

t = -AP"i, x = -pAF-^, p= -MF-^/2, (4.4) 

3 

where M = ^ fJ'i{t, x)p^ is a cubic polynomial in p with coefficients 

i=0 



^XJ 



/i3 = a{at — 26^) + ba^, ^2 = b{3at — 2bx) + ca^ — 2ac^ 

/ii = b{2bt - 3c^) + 2cat - act, fio = c{2bt - c^) - bct. 

"* 3 ~* 

Multiplying W by — P2, we obtain the field V 

t = A, i = pA, p = M/2. (4.5) 

For any point g* = (t*, x*) E £WH and any p G M.P there exists a unique geodesic passing 
through g* with given tangential direction p. However if g* is parabolic, this is not the case. 
Indeed, for any tangential direction p G M.P such that M{q^:,p) 7^ there exists a unique 
trajectory of (14. 5 p passing through the point {q*,p), a vertical line, which projects onto the 
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single point q^ in the (t, x)-plane. Thus, geodesies outgoing from g* must have tangential 
directions p such that M{q^,p) = 0, i.e., their 1-graphs pass through a singular point {q*,p) 
of the field V. 

Let g* G A and consider the equation M{q^,,p) = with respect to p. We shall assume that 
in a neighborhood of g* the curve A is regular and a(g*) 7^ 0. Then the quadratic polynomial 
F{q*tP) = ap^ + 26p + c has a unique root po{q*) = "^^rH? that is, the isotropic direction|f| 
A simple substitution shows that po is a root of the cubic polynomial M(g*,p). Assume that 
the isotropic direction pq is not tangent to the curve A at g*, i.e, {aAt — bAri.)\g^ ^ 0. 

Under the assumptions above, the cubic polynomial M{q^,p) has one or three real prime 
rootso Define Wq = {q & A, p = po{q)} and W^. = {q & A, p = p±{q)} where p±{q*) are 
the non-isotropic roots of M{q^^p) = 0, if they exist. The union of the three curves Wq, W^ 
is the set of singular points of V and coincides with its center manifold W^. The function 
F vanishes on Wq while F 7^ on W^. Thus the fields (14. 4p and (14. 5 p are connected by 
relation (ll.ip . where f = F and ^ = |. Since the field W is obtained from an Euler-Lagrange 
equation, conditions (II. 2p and (II. 3p follow from the identity D^ = 0, which is valid for all 
points except for the hypersurface F = {F = 0}. From Theorem [1] it follows that F is an 
invariant hypersurface of V. Hence the isotropic curves are geodesic lines (of zero length) in 
the pseudo-Riemannian metric (14. 3p . By construction W^ C F. 

Let (g,Po) ^ ^0- Clearly, the spectrum of the hnearization of V at (g,Po) contains the 
eigenvalue Ai = A^ + pqA^ 7^ 0. By Theorem [HI in a neighborhood of (g,Po) there exists a 
function a : Wq — )■ M such that the spectrum of the linearization of V at all points sufficiently 
close to (g,po) is cr(2, 1, 0). Computing, we easily get a = At + poA^;. Hence the germ of V 
at (g, po) is smoothly orbitally equivalent to 

i = 2^ + ^{0r]\ r] = r], ( = 0. (4.6) 

The normal form (14. 6 p can be further simplified. 

Theorem 9 The germs of the vector field V given by formula (14. 5 p at the singular points 
ilyPo) ^ W^o ^^^ (?!P±) ^ W^± ore smoothly orbitally equivalent to 

e = 2e, v = v, c = o, (4.7) 

e = e, V = -V, C = ^V, (4.8) 

respectively. 

Proof. To establish normal form (14. 7p it is sufficient to prove that the coefficient V'(C) in 
the normal form (14. 6 p is identically equal to zero. Indeed, the field (14. 6 p has the invariant 
foliation {{^, r/, (^) : ^ = c}, and the restriction to each leaf is a node with exponenci equal to 
2. The eigenvalue of largest modulus corresponds to the eigenvector ^ and the eigenvalue 

of smallest modulus corresponds to ^. 

Given an arbitrary point (g*,Po) ^ W^o, consider the restriction of the field (14. 6 p to the 
invariant leaf {{C,,ri,Q '■ C = C*} passing through (g*,po)- Integrating the corresponding 
differential equation dC,/dri = 2C,/ri + ip{(^^:)ri, we get the single integral curve rj = and the 
family of integral curves 

^ = c?7^ + (y9((^^,)?7^ In Ir^l, c = const, (4.9) 

^The light cone at a parabolic point consists of a unique isotropic line. 

■* If 5 is a surface embedded in 3D Minkowski space, these cases correspond to positive or negative Gaussian curvature of S 
in the Euclidean metric dx^ + dy^ + dz^. 

^ The exponent of a node (or saddle) is defined to be the ratio of the eigenvalue of largest modulus of the linearization field 
to the smallest one. 
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with common tangential direction ^ at 0. 

In the case ip{(*) = all curves of the family fl4.9p are parabolas, in the case (p{C*) ¥" 
they are C^-smooth, but not C^-smooth at 0. On the other hand, the previous reasoning 
shows that the germ of V at {q*,po) has at least one C°°-smooth integral curve: the vertical 
line (parallel to the p-axis). Simple calculation shows that the direction ^ in the initial 

coordinates (t, x,p) corresponds to the direction -^ in the normal coordinates (^, r], (). Hence 
family (14. 9p contains at least one C°°-smooth integral curve. This implies that v^(C*) = 0. 

The second statement of the theorem (the normal form (I4.8P ) follows from Lemma [3] and 
Theorem [71 validity of the condition \l/i(0,0) 7^ can be proved by direct calculation (see 
Theorem 2 in [I3]). ^ 

4.2 Metric of Klein type 

A natural generalization of the Klein metric on the (t, x)-plane is 

, 9 a dx^ + 2/3 dxdt + •y dt"^ ^^ . . 

ds^ = ^ 1 , riGN, (4.10) 

where the numerator is a positive definite quadratic fornix with coefficients a, 13, 7 smoothly 
depending on t, x. We study locally the geodesies of metric (14.1 OP passing through a singular 
point, i.e., a point of the axis {t = 0}. It is not hard to prove that in appropriate local 
coordinates on the (t, x)-plane the germ of metric (I4.10p simplifies to the form 

2 a dx'^ + 'y dt'^ ^^ / , -, . x 

^^ = ^ ' ^ ^ ^' (4.11) 

with smooth positive coefficients a{t,x) and 7(t,x). 

The geodesies of metric (14.110 are extremals of the Euler-Lagrange equation (14. ip with 
L = yF/t"", where F = ap^ + 7 > and p = dx/dt. The corresponding vector field W reads 

t = a7t-"F-i, y = a-fpt-''F-l, p = --t'^'-^MF-^ , (4.12) 

3 
where M = ^ fJ'i{t, x)p^ is a cubic polynomial of p with coefficients 

i=0 

II3 = a{tat - 2na) , /i2 = t(a^7 - 207^.), fii = t(2aa - a-^t) - '^na'j, fio = -t'^y^. 

Multiplying W by Z"'*'"'^, where / = tg'^^^ and g = F'^/^/(«7) > 0, we obtain the field V 

i = t, x = pt, p=-M/{2a^). (4.13) 

Fields (I4.12P and (I4.13p are connected by relation (II. ip . where / = tg'^^^ and r = n + 1. 
Conditions (ll.2p and (II. 3p are satisfied, W being obtained from an Euler-Lagrange equation. 
Theorem [1] implies that F = {/ = 0} = {t = 0} is an invariant plane for V. The restriction 
of the field V^ to F is parallel to the p-axis. Hence geodesies outgoing from a point g* = (0, x*) 
must have tangential directions corresponding to p such that M{q^,p) = 0, i.e., their 1-graphs 
pass through singular points of V. 

Given a point g* = (0,x^,), consider the equation M{q^,p) = with respect to p. Since 
M{.q*^p) = —'2nap{ap'^ + 'y), the cubic polynomial M{q^,p) has the only real root p = 0. The 



^ The case when the numerator is an indefinite (and non-degenerate) quadratic form was also studied 1141 . but for our 
present purposes it is sufficient to consider the positive definite case. 
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spectrum of the linearization of V at (g*, 0) is {n, 1, 0), and the x-axis is the center manifold 
{W^). Clearly, W^ C F and from Theorem [8] we get the following result. 

Theorem 10 The germ of the vector field (I4.13P at the singular point (g*,0) is smoothly 
orbitally equivalent to 

e = < + <^(C)r7", v = V, C = 0. (4.14) 

Unlike the case of geodesies in pseudo-Riemannian metrics, here the coefficient <f{C) is 
not necessarily zero. For instance, in the case n = 1 the condition v?(C*) = is equivalent to 
7x(0,x*) = 0, where C* corresponds to the point g* = (0,x*). Clearly, if ip{0) ^ then ^{Q 
simplifies to 1, if v^(C) has a finite order s at the origin then V2(C) simplifies to CJ^ . 

Example 6. Consider the Klein metric, given by formula (14.111) with a = 7 = 1 and 
n = 1. In this case vector field V given by (I4.13P has the normal form (14. 141) with V5(C) = 0, 
since 7j;(0,a;) = 0. Hence the restriction of the field V on each invariant leaf (given by the 
formula C = c in the normal coordinates) is a bicritical node. Thus the integral curves of V 
are C°°-smooth, and for each singular point q^, = (0, a;*) there exists a family of geodesies with 
common tangential directions p = and various 2-jets. Indeed, geodesies of the Klein metric 
passing through the point g* G A (here A is the absolute) are the circles (x — x,,,)^ + 1^ = i?^ 
and the straight line x = x*. 

4.3 Almost-Riemannian metric 

Let V^i, V2 be smooth vector fields on the (x, |/)-plane. Assume them to be collinear at the 
points of a regular curve A. The metric ds"^ having (Vi, V2) as orthonormal frame is well- 
defined, smooth and positive definite on the whole plane except for the points of A. Our 
aim is to study geodesies of the metric ds'^ passing through a point of the curve A. 

According to [T], for a generic pair Vi, V2 in a neighborhood of almost all points of A there 
exist local coordinates such that 

Vi = g^, V2 = 2xv-\x,y)—, (4.15) 

where v{x,y) is a smooth non-vanishing function. Points in A at which such a coordinate 
system does not exist form a discrete subset of A and will not be considered in the following. 
The required metric for fields (I4.15P is 

2 , 2 , ^^ . 2 v^ dy^ + djx^? 
ds =dx +—dy = . 

Substituting t = x^ and multiplying by the unessential constant factor 4, we get 

tP dxi^ -\- dt^ 
ds^ = '— , v{t,y) = vi{t,y) + Viv2{t,y), (4.16) 

where wi,2 are smooth functions defined by decomposition of the function v{x,y) into even 
and odd parts: v{x,y) = Vi{x'^,y) +xv2{x'^,y). Geodesies of the metric (I4.16P are extremals 
of Euler-Lagrange equation with Lagrangian L = y^F/t, where F = v^p^ + 1 and p = dy/dt. 
The corresponding vector field W in the (t, j/,p)-space reads 

i = ^t~^F-2, y = ^pt-^F-^, p=-t-^F~2M, (4.17) 
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__- 3 ^ 

where M = Y2 f^i{tjy)p^ is a cubic polynomial in p with coefficients 

Jl^ = v'^ — 2tv^Vt, /i2 = —2tVy, Jii = v — Atvt, yUo = 0. (4-18) 

Multiplying W hj f^, where r = ^, f = tg and g = {2/v)^F ^ 0, we obtain the field V 

i = 2vt, y = 2vtp, p = M. (4.19) 

Fields fl4.17p and fl4.19p are connected by relation ( II. ip . where the function f = tg is regular 
and r = ^. Nevertheless, in general we cannot apply Theorems [H |5] and [HI since the field V 

is not even C^-smooth. Indeed, the components of the field V depend on the function v{t, y) 
and its first-order derivatives, which are smooth only ii v{x, y) is an even function of x (see 
formula ( KWf ). 

Example 7. Consider the Clairaut-Liouville metric. This is an example in which the 
vector field (I4.19P turns out to be smooth, the function v{x, y) being even in x. For instance, 
in |3j the authors deal with the metric 



ds = dx H — dy 



g{x'^, y) , 2 2;2 dx'^ + fi'(3^^ v) dy'^ 

^dy = -2 



x^ x" 

where g is a, positive smooth function (x and y are standard angular coordinates on the 
sphere, the curve A = {x = 0} is the equator) After the change of variables t = x^ we get 
the metric (I4.16P with v = 2^/g{t/y), which leads to the smooth field (14.191) . 

To overcome the problem, we make the change of variable x^ = t in (I4.19p . This yields to 

X = XV, y = 2x'^vp, p = M, (4.20) 

3 

where M = J2 l^ii.^, y) p* and /ij(x, y) = fj.i{x^, y). The coefficients /ij are polynomials of the 

^ i=0 ^ 

function v{t,y) and its first-order derivatives (see formulas (I4.18P ). Note also that Vt appear 

in (I4.18P with the factor t, whence after the substitution x"^ = t the expression tvt becomes 

a smooth function of x, y. 

The first two components of the field (I4.20p vanish at x = 0. Given a point q = {0,y) 

consider the cubic equation M{q,p) = with respect to p. It reads v{q)p{{y{q)pY + 1) = 0. 

This equation has a unique real root po = 0. Recalling that p = dy/dt, the root Po = defines 

the unique admissible direction for geodesies passing through the point (0, y) on the (t, y)- 

plane. The corresponding direction on the (x, y)-plane is given by the relation dy/dx = 2xp 

which is also equal to zero. The spectrum of the linearization of the germ (I4.20p at {q,po) is 

(Ai, A2,0), where Ai = v{q) and A2 = Mp{q,po) = v{q). 

Theorem 11 The germ of the vector field (I4.20p at the singular point {q,po) is smoothly 
orbitally equivalent to 

e = e, V = V, C = 0. (4.21) 

Proof. By Theorem [6l the germ of the vector field (I4.20p at {q,po) is smoothly orbitally 
equivalent to normal form (I3.13P with X{() = 1 and / = 1. To establish normal form (I4.2ip 
it is sufficient to prove that the coefficient <f{C) in (I3.13P is identically equal to zero. 

Let A be the linearization of the vector field (I4.20p at the singular point {q,Po)- Consider 
the matrix A — XI, where A = v{q) is the double eigenvalue of A. Clearly, the value rg(A — A/) 



'^ In the case g{x^ ,y) = 1 this formula gives the well-known Grushin metric. 
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equals either 1 or 2 and it is an invariant of the field. Hence y^{0) = if rg(A — XI) = 1 and 
(f{0) 7^ if rg(A — A/) = 2. On the other hand, a simple calculation shows that rg(A — A/) = 1 
if and only if Mx{q,po) = 0. Recalling that M{q,p) = v{q)p{{v{q)pY + 1) and po = we get 
Mj;(g,po) = 0. This completes the proof, q 

From the normal form f l4.2ip it follows that vector field fl4.20p has an invariant foliation 
(given by C = const in the normal coordinates) such that each leaf intersects the center 
manifold W^ at a unique point. Hence, geodesies passing through the point g* = (0,?/*) 
on the {x, ?/)-plane are projections of integral curves lying in the corresponding leaf. The 
restriction of vector field f l4.20p to the leaf is a bicritical node, hence there is a one-parameter 
family of integral curves passing through the point g* = (0, y^,). This gives a family of smooth 
geodesies passing through the point g* with common tangential direction which coincides 
with Vi and % at the point q^. Moreover, the geodesies have the same 2-jet and different 
3-jets at g*. 

Example 8. Geodesies in the Grushin metric (which corresponds to the vector fields 
(I4.15P with v{x,y) = 2) have the form y{x) = y* + c~^ arcsin(cx) — c~^x\/l — c^x^, where c 
is an arbitrary constant. 
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